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ABSTRACT 



In a previous paper, we argued that the inversion of Stokes profiles applied to spectropolarimetric observations of the 
solar internetwork yield unrealistically large values of the inclination of the magnetic field vector (7). This is because 
photon noise in Stokes Q and U are interpreted by the inversion code as valid signals, that leads to an overestimation of 
the transverse component B±, thus the inclination 7. However, our study was based on the analysis of linear polarization 
signals that featured only uncorrelated noise. In this paper, we develop this idea further and study this effect in Stokes 
Q and U profiles that also show correlated noise. In addition, we extend our study to the three components of the 
magnetic field vector, as well as the magnetic filling factor a. With this, we confirm the tendency to overestimate 7 when 
inverting linear polarization profiles that, although non-zero, are still below the noise level. We also establish that the 
overestimation occurs mainly for magnetic fields that are nearly vertical 7 < 20° . This indicates that a reliable inference 
of the inclination of the magnetic field vector cannot be achieved by analyzing only Stokes I and V. In addition, when 
inverting Stokes Q and U profiles below the noise, the inversion code retrieves a randomly uniform distribution of the 
azimuth of the magnetic field vector (f>. To avoid these problems, we propose only inverting Stokes profiles for which 
the linear polarization signals are sufficiently above the noise level. However, this approach is also biased because, in 
spite of allowing for a very accurate retrieval of the magnetic field vector from the selected Stokes profiles, it selects 
only profiles arising from highly inclined magnetic fields. 

Key words. Magnetic fields - Sun: photosphere - Sun: surface magnetism - Stars: magnetic field 
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1. Introduction 

Spectropolarimetry, which is the study of the polarization 
properties of the light observed in spectral lines, is the 
most developed and widely used tool for retrieving the 
magnetic properties of the solar plasma and other astro- 
physical objects (Stenflo 2002; Mathys 2002 and references 
therein). The inference of the magnetic field vector from 
spectropolarimetric observations is performed through the 
inversion of the radiative transfer equation for polarized 
light (del Toro Iniesta 2003a; Bellot Rubio 2004; Ruiz 
Cobo 2007). However, as the module of the magnetic field 
vector decreases, the observed signals disappear below 
the level of the photon noise, making the inference of the 
magnetic field very difficult and plagued with problems and 
uncertainties. This is, for instance, the case for the areas 
on the solar surface referred to as the solar internetwork. 
Here, the polarization signals are so weak that for a long 
time the internetwork was thought to be void of magnetic 
fields. 

With the advancements in the sensitivity of the 
polarimeters and the higher spatial resolution of the 
observations achievable with adaptive-optic systems and 
large-aperture telescopes, polarization signals are now 
routinely detected everywhere in the internetwork (Lites 
et al. 1996). This demonstrates that these regions are truly 



pervaded by magnetic fields. Unfortunately, the signals 
are barely above the noise level, thus it has not yet been 
possible to fully characterize the magnetic field vector in 
these regions. This has led to a long-standing controversy 
about the distribution of the module of the magnetic field 
vector in the internetwork (Dommguez Cerdefia et al. 
2003; Socas-Navarro & Lites 2004; Sanchez Almeida 2005; 
Martinez Gonzalez et al. 2006; Asensio Ramos et al. 2007; 
Lopez Ariste et al. 2007; Socas-Navarro et al. 2008; see 
also references therein). 



More recently, further discrepancies have emerged 
about the angular distribution of the magnetic field vector 
(Lites et al. 2007, 2008; Orozco Suarez et al. 2007a, 
2007b; Martinez Gonzalez et al. 2008; Asensio Ramos 
2009; Stenflo 2010; Ishikawa & Tsuneta 2011; Borrero & 
Kobel 2011). In Borrero & Kobel (2011; hereafter referred 
to as paper I), we argued that the photon noise leads 
to a systematic overestimation of the inclination of the 
magnetic field vector. We reached this conclusion because 
when employing only vertical magnetic fields to synthesize 
Stokes profiles, these could be retrieved as highly inclined 
ones due to the sole effect of the photon noise. Since 
those tests were carried out with vertical magnetic fields, 
the linear polarization signals (Stokes Q and U) were 
zero. This limits the validity of the tests, as the inversion 
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code analyzes only linear polarization profiles featuring 
uncorrelated noise. It is therefore worth considering the 
case where small linear-polarization signals (due to non- 
vertical magnetic fields) can be hidden below the noise. 
This introduces correlations that could be employed by the 
inversion algorithms to retrieve useful information about 
the inclination of the magnetic field vector. This paper 
is devoted to addressing this particular situation. To this 
end, we employ uniform distributions of the magnetic field 
vector (Section 0) to produce a large database of Stokes 
profiles. Therefore, unlike paper I, the resulting Stokes 
profiles will have non-zero linear polarization profiles. 
After adding noise to these profiles, we apply an inversion 
algorithm (Section[3]) that attempts to retrieve the original 
distribution of the magnetic field vector. Our results are 
described in Sections 0] and [SJ Section [5] summarizes our 
findings. 



2. Synthesis with uniform distributions 

To study the effect of the photon noise and selection crite- 
rion, we performed several numerical experiments employ- 
ing a probability distribution function that is parametrized 
by the physical parameters of a Milnc-Eddington atmo- 
sphere (ME). These physical parameters are denoted as X 



X= [B,a,y los ,T] , 



(1) 



where B refers to the magnetic field vector, a is the 
so-called magnetic filling factor and denotes the fractional 
area within the pixel that is occupied by magnetized 
plasma, V\ os is the line-of-sight component of the velocity 
vector, and, finally, the vector T considers all thermo- 
dynamic quantities (source function, Doppler width, and 
so forth). More details about these parameters and the 
Milne-Eddington approximation can be found in Landi 
Degl'Innocenti (1992) and del Toro Iniesta (2003b). Here 
we write the probability distribution function for the M-E 
parameters X in the form 



P(X)dX = V 1 (B,a)V 2 (T, Mos)dBdadTdMc 



(2) 



where the thermodynamic and kinematic parameters, T 
and Vios, are assumed to be statistically independent of 
the magnetic ones B and a. The distribution function of 
the former parameters T-^T, Vi os ) is the one obtained from 
the inversion of map A in Borrero & Kobel (2011), so that 
we employ values that are representative of the quiet Sun. 
The term that contains the properties of the magnetic 
field, Pi(B,a), has been modeled as a slightly-modified 
uniform distribution, expressed as 



V 1 (B,a)dBda = 



1 



2tt 2 B, 



■R(B - B )dBd-fd(j)da 



(3) 



which gives the probability of finding a magnetic field 
vector B whose module is between B and B + dB, whose 
inclination (with respect to the observer's line-of-sight) is 
between 7 and j + dj, and whose azimuth (in the plane 
perpendicular to the observer's line-of-sight) is between 



(j> and <p + dip. On the right-hand-side of this equation, 
the term K(B — Bq) also refers to the complementary 
Heaviside function 



H(£ - B ) 



1, if .8 < Bo 
0, HB>B 



(4) 



where Bo is taken as Bq — 200 G in order to make our 
experiment representative of weak field regions in the solar 
surface. In addition to the value of Bq, we do not make 
any further attempts to employ a more realistic model for 
the solar internetwork since our aim is not to investigate 
particular distribution functions, but rather to study the 
effect of the inversion, photon noise and selection criteria 
in the inference of the magnetic field vector in a general way. 

Once we have the probability distribution function 
given by Eqs. [2][3]|H we solve the radiative transfer equa- 
tion in order to create a large database (= 2 x 10 6 ) of 
synthetic/theoretical Stokes profiles of the Fel 6302.5 A 
(ffeff = 2.5) spectral line. The effective Lande factor of the 
atomic transition calculated under the LS coupling scheme 
is indicated by g e g. The Stokes profiles are synthesized 
with the VFISV (Very Fast Inversion of the Stokes Vector) 
inversion code of Borrero et al. (2010). To these profiles, 
photon noise is added as a normally distributed random 
variable with a standard deviation of a = 10~ 3 ,3 x 10 -4 . 
These two values mimic the noise levels of maps A and 
C in Borrero & Kobel (2011). Once the noise is added, 
the resulting Stokes profiles are taken as real observations 
and inverted in order to retrieve the original atmospheric 
parameters X. This step is described in the next section. 



3. Inversions of synthetic Stokes profiles 

The Stokes profiles for Fel 6302.5 A synthesized in the 
previous section are now inverted with the same VFISV 
code, but now running in inversion mode instead of synthe- 
sis mode. At the start, VFISV solves the radiative transfer 
equation for polarized light in the Milne-Eddington (ME) 
approximation using a set of initial values for the physical 
parameters: Xo. The solution of the radiative transfer 
equation yields theoretical or synthetic Stokes profiles 
that are compared to the ones synthesized in the previous 
section (those where photon noise had been added). The 
initial values of Xo are then iteratively modified until 
the best possible fit between the theoretical/synthetic 
and observed Stokes vector is reached. The final Xf that 
achieves the best fit is then assumed to be the real one 
present in the solar atmosphere. In this work we will 
focus mainly on the magnetic field vector B and magnetic 
filling factor a. For more details about how the inversion 
is performed, additional free parameters of the inversion, 
and treatment of the filling factor, we refer the reader to 
Borrero et al. (2010) and Borrero & Kobel (2011; paper 
I). We note that, since we are using the same type of 
atmospheric model in the synthesis and inversion (ME 
atmospheres), the experiments carried out in this paper do 
not address the systematic errors introduced by the choice 
of the wrong atmospheric model, such as those introduced 
when inverting asymmetric Stokes profiles using a ME 
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atmospheric modeQ- 

Once B, 7, <p and a had been retrieved from each 
Stokes profiles of the original database, we applied two 
different selection criteria. In the first criteria, we select for 
the analysis profiles in the database where the maximum 
of the absolute value (for all wavelengths) in any of the 
polarization signals (Stokes Q, U, or V) is larger than 4.5 
times the noise level max|<5(A), U{X), V(X)\ > 4.5er. This 
is equivalent to selecting pixels where the signal-to-noise 
ratio (S/R) in the polarization profiles is 4.5 or better, 
i.e., S/R > 4.5. Hereafter we refer to this criteria as 
S/R q uv In the second criteria we select those profiles 
within the database where the maximum of the absolute 
value (for all wavelengths) in the linear polarization signals 
(Stokes Q, U) is larger than 4.5 times the noise level 
max|Q(A), U(X)\ > 4.5a. This criteria will be referred to 
as S I -Rqu • 

The reason behind the choice of these two different 
selection criteria is the following. The first criterion (selec- 
tion of profiles where S/R quv > 4.5) was adopted so we can 
draw parallelisms with both the results presented in paper 
I and Orozco Suarez et al. (2007a, 2007b). Owing to the 
intrinsic amplitude of Stokes Q, U, and V, this criterion 
selects mostly pixels where only the circular polarization 
V possesses a S/R > 4.5. The actual percentages depend 
on the distribution of the magnetic field vector and the 
level of noise. For instance, taking a uniform distribution 
of B, 7, and <j) (Equation and considering a noise-level 
of g = 10~ 3 , only 21.2 % of the profiles selected with the 
Sj -Rquv-criterion possess linear polarization signals (Stokes 
Q and U) with peak-values above 4.5er. In the remaining 
78.8 % of the profiles, only Stokes V is above 4.5<r. The 
21.2 % increases up to 39.1 % when considering a noise 
level of a = 3 x 10~ 4 . In paper I, the maps that had 
equivalent levels of noise featured different percentages of 
profiles: in the map with a = 10 ~ 3 , only about 8.8 % of 
the profiles selected with the iSyi?q U v-criterion had linear 
polarization profiles above 4.5cr. This number increased 
to 38.8 % in the map with a = 3 x 10~ 4 . As already 
demonstrated in that paper (see also Martinez-Gonzalez 
et al. 2011), the inversion of profiles selected with the 
S/ iiq UV -criterion yields values for the inclination of the 
magnetic field vector, 7, that are largely overestimated. 
This happens because most of the Stokes profiles have 
linear polarization signals that are either below or close 
to the noise level. This effect decreases as the noise in 
the polarization profiles goes down so that the linear 
polarization signal rises above the noise. However, since 
both Stokes Q and U remain below 10 -4 for transverse 
fields up to B± w 20 — 30 G, the overestimation in the in- 
clination remains significant even at such low levels of noise. 

In paper I, we proposed, as an alternative solution, 
selecting for the inversion only pixels where the peak in the 
Q or U signals is at least 4.5 times higher than the noise 



1 Other sources of systematic errors that are not being consid- 
ered, as they are beyond the scope of our study, are the effects 
of the spatial resolution, spectral resolution, etcetera. For in- 
stance, Borrero et al. (2007) carried out such a study for the 
particular case of the Helioseismic and Magnetic Imager (HMI) 
instrument. 



i.e., where the S/R is higher than 4.5 in linear polarization. 
We also followed this approach in applying the second 
criterion (see also Asensio Ramos 2009). In paper I, 
we anticipated that this approach has the advantage of 
retrieving more trustworthy values of 7, B±, and B but we 
did not quantify it. 

With the results from the selected profiles using the 
aforementioned selection criteria, we constructed his- 
tograms for the magnetic field strength B, the inclination 
of the magnetic field vector with respect to the observer's 
line-of-sight 7, the projection of the magnetic field vector 
along the observer's line-of-sight B\\ = -BCOS7, the projec- 
tion of the magnetic field vector along the perpendicular 
direction to the observer's line-of-sight: B± = Bsmj; 
azimuthal angle of the magnetic field vector in the plane 
perpendicular to the observer: 0, and finally the magnetic 
filling factor a. Results for a = 10 -3 and S/R quv are 
displayed in Figure [TJ whereas Figure [5] presents the 
results for a = 3 x 10~ 4 and S/R quv . Likewise, Figures [3] 
and |4] present the results for these same levels of noise, 
respectively, but employing the S/R qu selection criterion. 
The panels in each of these four figures correspond to the 
following histograms: a) the component of the magnetic 
field vector along the observer's line-of-sight Bn; b) the 
component of the magnetic field vector perpendicular the 
observer's line-of-sight B±; c) the module of the magnetic 
field vector B; d) the inclination of the magnetic field 
vector with respect to the observer's line-of-sight 7; e) the 
azimuthal angle of the magnetic field vector in the plane 
perpendicular to the observer's line-of-sight (f>; and finally 
f) -panels for the magnetic filling factor a. In all panels, 
solid-black lines represent the initial distribution employed 
in the synthesis of the Stokes profiles (see Eqs. HE]) 
including all 2 x 10 6 Stokes profiles of the distribution. 
As required, the distributions for B, 7, </>, and a are 
uniform (equal probabilities). We note that this is not 
the case for B\\ = Bcos'f and B± = Bsinj because the 
cosine and sine functions introduce non-uniformities in 
the distribution. In addition, dashed-black lines represent 
the original distribution of the physical parameters but 
considering only the profiles that are selected with corre- 
sponding selection criteria, namely S/R quv in Figs. [TJH 
and S/R qu in Figs. [SB Finally, the solid-red lines present 
the distributions inferred from the inversion of the Stokes 
profiles but employing the same Stokes profiles that were 
used to construct the dashed-black histograms. In the 
following we comment on some general features that will 
assist us in interpreting our results: 



— The area under the solid-black curve is always one, 
while the area under the dashed-black and solid-red 
curves are both normalized to the quotient of the 
number of profiles selected with each particular criteria 
and the total number of profiles synthesized (2 x 10 6 ). 
Therefore, one can interpret the difference between 
the solid-black and dashed-black curves as the loss 
of information caused by selection criterion, whereas 
the differences between dashed-black and solid-red are 
caused by the inversion algorithm. 

— If the solid-red curve lies below the dash-black one 
for a given interval of the physical parameter (e.g. 
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A7), it means that the number of Stokes profiles that 
were synthesized employing that range of values (e.g. 
A7) are underestimated by the inversion. Likewise, 
whenever the solid-red curve lies above the dash-black 
one for a given interval, the inversion overestimates 
the number of Stokes profiles that were synthesized 
employing the range of values of the physical parameter 
given by the same interval. 

— The area under the dashed-black lines in each panel 
increases as the noise level decreases. This can be 
realized by comparing Figure [I] (a = 1CP 3 ) with 
Figure [5] (a = 3 x 10~ 4 ), and also comparing Figure 
[3] ( a = 10- 3 ) with Figure H (a = 3 x I(T 4 ). This 
happens as a consequence of a larger number of Stokes 
profiles fulfilling the requirement imposed by a given 
selection criterion when the noise is reduced. Since 
the S/R qu is more stringent than S/R quv (because 
linear polarization signals are usually much weaker 
than the circular polarization signals), the area under 
the dashed-black curves is of course much larger in 
Figs. [Q[2] than Figs. EH 



4. Discussion 

We now discuss qualitatively the effects of photon noise 
and selection criteria for each physical parameter indi- 
vidually. We first introduce the approximate dependences 
of the Stokes profiles on the three spherical coordinates 
of the magnetic field vector (B, 7, <j>) and the magnetic 
filling factor a, in the weak field approximation (Landi 
Degl'Innocenti 1992) 

V oc aB cos 7 

Q oc aB 2 sin 2 7 cos 2<j) (5) 
U oc aB 2 sin 2 7 sin 2<j> 

4.1. Line-of-sight component of the magnetic field: B» 

As mentioned in the previous Section, the Sy-Rq UV -criterion 
selects a portion of the initial Stokes profiles (dashed-black 
lines in Figs. [T^,-[2ji) that more closely resembles the 
original B\\ distribution (solid-black lines) than profiles 
selected by the S , /i? qu -criterion (dashed-black lines in 
Figs. [3hr|lk). This is because the S'/iiquv-criterion selects a 
much larger sample of Stokes profiles than S/R qu , hence 
the distribution obtained with the former criterion must 
more closely resemble the original distribution employed 
to construct the database of Stokes profiles (solid-black 
lines). 

We also note that the distribution of the component of 
the magnetic field vector along the observer's line-of-sight 
.By is very accurately retrieved by the inversion in all 
instances, regardless of the selection criterion and either of 
the noise levels (see Figs.QJi-rjJi). This conclusion is reached 
from the fact that solid-red lines coincide almost perfectly 
with the dashed-black lines in all figures. In addition, the 
solid-red curves shift towards lower values of B\\ as the 
noise decreases, which indicates that lower levels of noise 



allow us to correctly infer smaller values of the line-of-sight 
component of the magnetic field. 



4.2. Transverse component of the magnetic field: B± 

For the same reason given above in the case of Bn, the 
SV-Rquv-criterion selects a portion of the initial Stokes 
profiles (dashed-black lines in Figs. [TJd-I^p) that resembles 
much more closely the original B±_ distribution (solid- 
black lines) than profiles selected by the 5/i? qu -criterion 
(dashed-black lines in Figs. EJd-HJd). In addition, the distri- 
bution obtained with the S/Rqu-criterion mostly provides 
information about the large values of the transverse 
component of the magnetic field (B± « 100 — 150 G), 
whereas the S'/iiquv-criterion also selects profiles arising 
from weak transverse magnetic fields (B± < 50 G). The 
reason for this is that much larger values of B± are needed 
to produce linear polarization profiles above the 4.5cr-level. 

The inversion code very reliably retrieves the distri- 
bution of B± employing the »S/Rq U -criterion: we note 
the almost perfect match between the solid-red curves 
and dashed-black ones in Figs. [3Jd-|Ud. However, this is 
clearly not the case for the S , /i? quv -criterion, since here the 
inversion underestimates the original distribution in the 
region where B± < 25 — 40 G (depending on the noise), 
but overestimates it in the region above this threshold. It 
is also noteworthy that the results from the inversion of 
the Stokes profiles selected with the 5/ i? quv -criterion show 
a Maxwellian-like distribution, where the peak appears 
progressively at lower values of B± as the noise decreases: 
fiP eak « 60 G (Fig. [Do; a = IO- 3 ) and B^ w 35 G 
(Fig. ^jp; a — 3 x 10~ 4 ). This is a consequence of having a 
higher sensitivity for lower levels of noise. 



4.3. Total magnetic field strength: B 

Since B = \J B 2 + B\ , and since the inversion of Stokes 

profiles retrieves very reliable distributions for B«, any 
mismatch between the original distributions of B (dashed- 
black) and the inferred ones (solid-blue) in Figs. [1]>|3J:, can 
be attributed to the same issues as the inference of B± 
(see above). In particular, as for B±, the 5/i? qu -criteria 
allows us to retrieve the distribution of the module of the 
magnetic field vector B reliably, while the S/i? quv -criteria 
underestimates the real distribution for values B < 35 — 70 
G (depending on the noise) but overestimates it above this 
threshold (solid-blue lines in Figs.[TJ: and[5£). 



4.4. Inclination of the magnetic field vector: 7 

As happens with Bu and B±, the Syi?q UV -criterion selects a 
set of Stokes profiles (dashed-black lines in Figs[TJi-[2l) that 
are more representative of the initial uniform distribution 
in 7 (solid-black lines) than the Stokes profiles chosen 
with the 5/i? qu -criteria (dashed-black lines in Figs EJi-Hti) - 
The Sy-Rq UV -criterion selects all sorts of inclinations, with 
a slight preference for more longitudinal magnetic fields. 
The Syi?q U -criterion selects however only highly inclined 
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fields (7 — > 90°). This happens because for weak magnetic 
fields, the only way to have linear polarization signals 
above the 4.5u-level is for these magnetic fields to be 
highly inclined. This can also be understood in terms of 
the kind of distributions obtained in Bu and in B± (panels 
a) and b), respectively) by each selection criteria and 
remembering that 7 = tan _1 (-B_i_/.B||). 

By comparing the dashed-black and solid-red lines 
in Figs. [Ui-|3Ji, we can conclude that the inversions of 
the Stokes profiles selected with the S , /i? quv -criterion 
yields a distribution of 7 that is underestimated for more 
longitudinal magnetic fields, but overestimated for more 
transverse ones. We note that the turning point between 
the underestimation of vertical magnetic fields and the 
overestimation of horizontal ones occurs at lower values of 
7 as the noise decreases: 7 ~ 22° for a — 10 -3 but only 
at 7 ~ 17° for tr = 3x 10~ 4 . Interestingly, the probability 
distribution function around 7 = 90° (where magnetic 
fields are perpendicular to the line-of-sight) is reliably 
inferred with this criterion. Altogether, these findings 
confirm our previous results (see paper I), for which 
we concluded that the inversion of the Stokes profiles 
obtained with the «Syi? quv -criterion will correctly infer 
very inclined magnetic fields, whenever these are present, 
but unfortunately also interpret as very inclined magnetic 
fields those that are mostly aligned with the observer's 
line-of-sight. In contrast, the inversion of pixels selected 
with the Sy-Rqu-criterion retrieves almost perfectly the 
distribution of the selected pixels, as can be seen by com- 
paring the dashed-black and solid- red curves in Figs.[5Ji-|3Ji. 



4.5. Azimuth of the magnetic field vector: (f> 

Here, both the S/R quv and S , /i? qu -criteria select a set 
of Stokes profiles that are representative of the original 
uniform azimuthal distribution of the magnetic field 
vector, 0, as the selected distributions are also very close 
to uniform (compare the dashed-black and solid-black lines 
in Figures [T^-|4^) . 

Remarkably, the inversion of the selected Stokes profiles 
for both values of the photon noise, er, and both selection 
criteria is able to retrieve the original distribution of the 
selected profiles, as the solid-red lines in Figures [lb-H^ 
match the dashed-black ones. To a first approximation 
(Auer et al. 1977; Jeferries & Mickey 1991), the az- 
imuthal angle of the magnetic field vector is given by 
4> = (1/2) tan _1 ([//Q). Thus, that the inversion of the 
Stokes profiles selected with the 5/-R quv -criterion also 
retrieves the correct distribution for <f> comes as a rather 
surprising result, since here most of the selected Stokes Q 
and U profiles are below the 4.5cr-level (see Section [3]). We 
address this particular point in Section [51 



4.6. Magnetic filling factor: a 

In the case of the magnetic filling factor a, neither the 
S/Rquv nor the S/R qu criteria (dashed-black lines in 
Figs. DJ-HF) are able to recover correctly the original 
uniform distribution (solid-black lines). In particular, most 



of the Stokes profiles arising from low values of the filling 
factor, a < 0.4, are neglected by both selection criteria. 
This is because the polarization signals scale linearly 
with a (see Equation [5J and therefore, small values of 
the magnetic filling factor yield polarization signals that 
are below the threshold employed in the selection. The 
situation is aggravated in the case of S/R qu (Figs. [30), 
where even larger values of a are needed owing to Stokes 
Q and U being intrinsically smaller than Stokes V. 

As far as the inversion is concerned, it is clear that the 
distribution of the magnetic filling factor of the selected 
profiles is very well-retrieved (solid-red lines in Figs. QJ-JU) 
for both noise levels and when employing both selection 
criteria. Interestingly, as happened with 7 (although 
to a much smaller extent), the inversion code slightly 
overestimates the selected distribution (dashed-black) for 
a < 0.5, but underestimates it above this value. 

We have so far only discussed the ability of the inversion 
code to retrieve the correct distribution for the magnetic 
parameters (three components of the magnetic field vector 
and filling factor). Although this is clearly an important 
question, it does not provide much information about the 
reliability of the inversion in individual cases. To address 
this point, we display in circles in Figure [5j the mean value 
of the differences between the original physical parameters 
of the selected Stokes profiles and the inferred ones through 
the inversion X sel — X lnv . These are denoted as AX i7 with 
Xi being any of the components of X in Equation [TJ which 
we refer to as bias in Xi. In addition to this, we also plot 
the standard deviation around the mean in the inference of 
the physical parameters as the dashed-color lines, which 
we refer to as a x (not to be confused with the photon 
noise a). For a physical parameter to be well-constrained, 
it is important that the mean is centered around zero 
(otherwise a systematic bias occurs) and that the standard 
deviation is small. 

We note that the size of the intervals in which the 
bias and standard deviation are calculated is not constant. 
This happens as a consequence of the lack of profiles in 
certain ranges. For instance, the S/R qu criterion barely 
selects profiles arising from a magnetic field vector where 
B± < 50 G (see red and yellow circles in Figs. |3b-[lb). This 
makes it necessary to increase the interval of B± (around 
this range) so that a sufficient number of profiles can be 
employed to obtain a statistically meaningful AB± and 
<jb± (see green and yellow circles in Figure [5b). This is not 
needed as often in the S , /i? quv -criterion because it selects 
many more profiles than the S'/i? qu -criterion (blue and red 
circles in Fig. [3b). 

By inspection of Figure [5j we note that, for the same 
level of noise and with the exception of the magnetic 
filling factor a, the standard deviation in the retrieval 
of the physical parameters is always much smaller when 
employing the S'/i? qu -criterion than the S'/i? quv -criterion. 
By comparing the standard deviations with different 
levels of photon noise, we realize that for most physical 
parameters, the improvement achieved by a decrease in 
the photon noise (from a = 10~ 3 to 3 x 10~ 4 ) is smaller 
than the improvement achieved by using the S/R qu 
instead of the S'/iiquv-criterion. This, generally makes the 
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dashed-green curves (a = 10~ 3 and Syi?q U -criterion) lie 
below the red-dashed ones (a = 3 x 10~ 4 and S/R quv - 
criterion). Interestingly, this does not apply to all physical 
parameters: in the case of B\\ and a, the photon noise plays 
a more important role than the selection criteria itself, 
as the red curves lie below the green ones. The relative 
importances of the selection criteria and the photon noise 
logically depends on the values of the noise considered, 
thus we cannot conclude that in general the former is more 
important than the latter (or the other way around) for 
some particular physical parameter. 

We now discuss the errors in the case of the S/R quv - 
criterion. Here, the retrieval of B± (blue and red colors in 
Fig. [5}d) has a large systematic bias towards _B" 1V > B s f l , 
such as AB ± < -75 G, for Bf l < 50 G. Furthermore, 
in this same region, the standard deviation is as large 
as &b± ~ 20 — 50 G. These numbers decrease as B± l 
increases, such that AB± ps —15 G and <jb± ~ 10 G for 
B s f l > 150 G. As happened in the case of the distributions 
in Figs. HJ:-|lt: the bias and standard deviations in the 
module of the magnetic field vector B (Figure [5k ) mimic 
those of B± (Figure Eb) . In the case of Bu (blue and red 
colors in Fig. [SJi), both the bias and standard deviation 
are always quite small < 5 G, owing to the Stokes V 
profiles typically being much larger than Q and U for 
the same values of B» and B±. For the inclination of 
the magnetic field vector with respect to the observer's 
line-of-sight, 7 (blue and red colors in Fig. [5ji), the bias is 
systematic towards more inclined magnetic fields: A7 < 
if 7 scl < 90°, but A7 > if 7 s01 > 90°. The absolute value 
of the bias is indeed larger for magnetic fields aligned 
with the observer's line-of-sight, where IA7I ss 30 — 40° 
if 7 scl ps 0, 180°. Interestingly, the bias decreases as the 
magnetic field vector becomes more perpendicular to the 
line-of-sight, i.e., A7 = for 7 s01 ps 90°. The behavior of 07 
is very similar (at all values) to that of the absolute values 
of bias IA7I. That the bias and standard deviation decrease 
as 7 — > 90° indicates that when the magnetic field vector 
is not completely aligned with the observer's line-of-sight, 
the signal that appears in Q and U helps the inversion code 
to extract some information about the inclination of the 
magnetic field vector. This inference is however negatively 
affected by the photon noise systematically overestimating 
7. This agrees with our results in paper I, where we found 
that it was impossible to distinguish between a vertical 
and horizontal magnetic-field vector when Stokes Q and U 
are below the 4.5u-level. Figure [SJ3 displays A<j> and as 
a function of 7 scl instead of <f> . In this case, the bias in 
the determination of the azimuthal angle of the magnetic 
field vector is almost negligible for all possible values of 
7 scl , where A<fi ps 0°. However, the standard deviation is 
large (er^ ~ 40 — 50°) when the magnetic field vector is 
aligned with the observer's line-of-sight (7 — > 0, 180°), but 
decreases as the magnetic field vector becomes more and 
more inclined with respect to the observer's line-of-sight 
(7 — > 90°). This happens as a consequence of the linear 
polarization profiles vanishing when the magnetic field 
becomes aligned with the observer such that Q, U cx sin 2 7 
(see Eq. H2). 

Compared to the S , /i? quv -criterion, there is no bias in 
the retrieval of all parameters (except for the magnetic fill- 



ing factor) when employing the S , /i? qu -criterion AX; — > 0. 
In addition, the standard deviation in the retrieval of the 
physical parameters, ct x , greatly decreases when employing 
the 5/i? qu -criterion, as indicated by the green and yellow 
colors in Figure [5] With this criterion AB± ps G and 
<tb ± ~ 10 G, even for Bj_ < 50 G. The retrieval of B\\ also 
improves, (such that crgy < 3 G) from that achieved with 
the iSY-Rquv-criterion. Finally, this translates into smaller 
standard deviations for B, 7, and <f>: gb < 5 G, <r 7 < 2°, 
and cx (p < 2°. 



5. Can inversion codes retrieve the distribution of 
the orientation of the magnetic field vector from 
correlations within the noise ? 

5.1. Azimuth of the magnetic field vector: cf> 

In Section 14.51 we have seen that the retrieved distribution 
of the azimuthal angle of the magnetic field vector on 
the plane that is perpendicular to the observer's line-of- 
sight (0) matches excellently well the original uniform 
distribution for the selected pixels. In the case of the 
5/i? quv -criterion, this seems counter-intuitive since many 
of the selected Stokes profiles only have sufficient signal 
in Stokes V, while Q and U are usually dominated by 
photon noise. One possible explanation for this is that the 
inversion does not properly infer <j) but instead yields a 
random distribution of values owing to the lack of informa- 
tion in the linear polarization profiles. If this distribution 
happens to be uniformly distributed, then the inferred 
distribution of azimuths matches the original one, even 
if the inversion fails for each Stokes vector individually. 
The second explanation is that, given that the number 
of profiles with large enough signals in Stokes Q and U 
is non-negligible, the inversion indeed properly retrieves 
cj>. In particular, in Section [3] we mentioned that 21.2 % 
(a = 10- 3 ) and 39.2 % (a = 3 x 10~ 4 ) of the profiles 
selected with the Sy.Rq UV -criteroii, had linear polarization 
signals above the 4.5cr-level. A final possible explanation is 
that the inversion is able to retrieve the correct values of <f> 
even when Stokes Q and U are dominated by noise. This 
could happen if the hidden signal introduces a sufficiently 
strong correlation in the noise, such that the minimization 
process can properly retrieve the azimuth of the magnetic 
field vector. To investigate which of these three possibilities 
are responsible for our results in Section [4J we repeated 
the synthesis (Sect. [5]) and inversion experiments of the 
spectral line Fe I 6302.5 A, but employing the probability 
distribution function 



7>i(B, a)dBda = 6 (<t>--) H(B - B )dBdjd(l)da (6) 

ttBq V 8 / 

which is very similar to Equation [3] with the exception 
that the distribution of the azimuthal angle is now a 
(5-Dirac centered at 22.5° = tt/8 rad. This value is cho- 
sen such that Q oc sin 2 7 cos 20, and U cx sin 2 7 sin 20, 
have the same amplitudes (see Eq. The distribution 
given by Equation [6] is employed to produce synthetic 
Stokes profiles, which are then inverted after adding 
photon noise with levels of a = 10 -3 and 3 x 10~ 4 . The 
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resulting distributions, for both levels of noise and both se- 
lection criteria (S/R quv and S/R qu ), are shown in Figure[6] 

Figure [5] shows that, on the one hand, the S/R quv - 
criterion (solid lines) has an almost uniform background 
distribution. This is caused by the inversion of Stokes 
profiles whose linear polarization signals are so weak that 
they carry no information about the azimuth, thus the 
inversion retrieves random results that happen to be uni- 
formly distributed. We know this because this background 
uniform field disappears if we apply the Syi? qu -criterion 
(dashed lines), which selects only Stokes Q and U profiles 
that are 4.5 times above the noise level. This is therefore 
one of the reasons for the excellent results in Figures [TJ:-[2t. 

In Figure [6j the 57-Rquv-criterion also features (super- 
imposed on the uniform background distribution) a clear 
reminiscent signature of the original 5(<f> — 7r/8) distribution 
(indicated by the vertical black arrow). This indicates 
that the inversion code is able to retrieve, to some extent, 
information about the azimuthal angle of the magnetic 
field vector. However, we still do not know whether <j) is 
correctly retrieved in only 21.2/39.2 % (a = 10 -3 , 3 x 1CT 4 , 
respectively) of the Stokes profiles that have sufficient 
signal in the linear polarization or, whether the inversion 
can properly infer (f> even when the linear polarization 
signals are below the 4.5a-level. To answer this question, 
we devised a third selection criterion, S/R 



0.04 



no— qu ) 



for 

which we selected all profiles where the linear polarization 
signals were below the 4.5<7-level. The distribution of 
<f> obtained with this criterion (dotted-dashed lines in 
Figure shows an almost uniform distribution, with a far 
smaller peak at <j> = ir/8 than before. This indicates that 
when the signals are dominated by noise, the correlation 
hidden below the noise does not provide sufficient informa- 
tion about the azimuthal angle of the magnetic field vector. 



5.2. Inclination of the magnetic field vector: 7 

In paper I, we demonstrated that a magnetic field 
vector that is aligned with the observer's line-of-sight 
(7 = 0, 180°) can be easily confused with a more per- 
pendicular one (7 w 90°) owing to the effect of the 
photon noise. In that experiment, the signals in the linear 
polarization were solely due to the photon noise. We did 
not study whether the inversion code could retrieve useful 
information about 7 from the Stokes Q and U profiles 
that, in spite of being below the 4.5cr-level, introduce a 
correlation within the noise. To study this possibility, 
we carry out a similar experiment to the one we have 
performed in Section 15.11 but focusing instead on the 
inclination of the magnetic field vector with respect to the 
observer's line-of-sight 7. To this end, we prescribe a the- 
oretical distribution of the magnetic field vector in the form 



Pl(B, a)dBda = ^]f 6 ~ \ ) H ( S ~ S o )d J Bd 7 d0da (7) 

which is identical to the distribution in Equation [6] but 
where the roles of <\> and 7 are exchanged. That is, 
the distribution of the azimuth is now uniform and the 
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Fig. 6. Histograms of the inferred azimuthal angle of the 
magnetic field vector, <f>, when employing Equation|B]in the 
synthesis of Stokes profiles. The original distribution 5{4> — 
7r/8) is indicated by the vertical black arrow located at <j> = 
22.5°. The inferred histo grams using S / R qU v , S / R qu , and 
S/R no -q U criteria are indicated by the solid, dashed, and 
dotted-dashed lines, respectively. Blue curves are obtained 
with a photon noise a = 10~ 3 , whereas red ones correspond 
to a = 3 x 1CT 4 . 



distribution of the inclination is now a (5-Dirac function 
centered at 7 = 45° = 7r/4 rad, which is chosen to ensure 
that B\\ = B cos 7 and B± — Bsiw-f (see Eq. [6]) are equal. 
As we did in Section 15.21 we add to the resulting Stokes 
profiles synthesized with Equation [7] photon noise with 
the values a = 10 _3 ,3 x 10 -4 . We also apply the three 
selection criteria mentioned above of S/R quv , S/R qu , and 
S/R no - qu . The inferred histograms of 7 are displayed in 
Figure 

With the first selection criterion, S/R quv , we realize 
that the inversion code retrieves a distribution of 7 (solid 
lines in Fig. [7]) with a clear peak at the original value 
of 7 = 45°. This distribution has however an extended 
asymmetric tail towards higher values of 7. This indicates 
that there is a tendency to retrieve field that are more 
inclined than they originally are, as a consequence of 
inverting Stokes profiles where Q and U are not 4.5 
times above the noise level. This becomes clear when 
we consider the 57-Rqu-criterion (dashed lines in Fig. [7]), 
which ensures that the linear polarization profiles are 
above the 4.5<r-level, since in this case the tail towards 
larger values of 7 disappears. What are the properties of 
the peak at 7 = 45° obtained with the 5/i? quv -criterion 
? Is it a consequence of the inversion code being able to 
properly retrieve the inclination from the weak linear- 
polarization signals that are below the noise level but 
introduce a correlation into the noise ? Or does the 
peak contain, as in the case of <f> (Sect. 15. ip and since 
a large portion of the Stokes profiles were selected with 
the 5/i? quv -criterion, strong linear-polarization signals ? 
This can be answered by the 5 l /i? no -q U -criterion, which 
only selects Stokes profiles for which Stokes Q and U are 
below the 4.5tr-level (see dotted-dashed lines in Fig. [7J. 
This selection criteria shows an even more pronounced 
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Fig. 7. Histograms of the inferred inclination angle of the 
magnetic field vector, 7, when employing Equation [7J in 
the synthesis of Stokes profiles. The original distribution 
£(7— t/4) is indicated by the vertical black arrow located at 
<j) = 55°. The inferred histograms using S/R quv , S/R qu , and 
S f /i? no _ qu criteria are indicated by the solid, dashed, and 
dotted-dashed lines, respectively. Blue curves are obtained 
with a photon noise a = 10 -3 , whereas red ones correspond 
to a = 3 x 1CT 4 . 

tail towards larger values of 7 than the S/i? quv -criteria. 
This is a consequence of many more Stokes profiles 
containing no information about 7 owing to very noisy 
linear-polarization signals. However, the peak at 7 = 45° 
does not completely disappear, which indicates that, unlike 
the case of 0, the inversion code is able to partially use the 
correlation hidden in the noise as a means of determining 7. 

We note, however, that this does not mean that the 
inversion code is able to retrieve the correct distribution 
of 7 by employing the S , /i? no -q U -criterion (Stokes Q and 
U profiles below the 4.5<7-level) when the distribution is 
a more general one than the <5-Dirac function employed 
in Equation [7] To emphasize this statement, it suffices to 
point out that not even the S , /i? quv -criteria, where 20-40 
% of the Stokes profiles have linear polarization signals 
above the 4.5tr-level, is able to do so (see Figs. QJi-[2Ji)- 
As demonstrated in Section [4J the only way of retrieving 
the original distributions of the magnetic field vector is to 
ensure that the linear polarization is above the 4.5<7-level, 
that is, to apply the S/i? qu -criterion. 



6. Conclusions 

Wc have carried out several numerical experiments in which 
a large number of Stokes profiles have been synthesized 
using uniform distributions of the three components of the 
magnetic field vector B, 7, and <f>, and the magnetic filling 
factor a (Equation [3]) . To these Stokes profiles, we have 
added photon noise to two different levels of a = 10~ 3 , 
and a = 3 x 10~ 4 . We have then inverted these profiles in 
order to retrieve the original distributions employed in the 
synthesis. This has been done in two different ways. In the 
first one, we have selected the pixels where the S/R ratio 



in any of the polarization signals (circular or linear) are 
equal to or larger than 4.5, the so-calld S , /i? quv -criterion. 
In the second case, we have selected the pixels where the 
S/R ratio in the linear polarization signals is equal to or 
larger than 4.5, which is the 5/i? qu -criterion. The former 
criterion, S/R quv , selects Stokes profiles that arise from 
all possible ranges of B, 7, and <j> (dashed-black curves 
in Fig. OH]), whereas the 5'/i? qu -criterion selects Stokes 
profiles that arise from highly inclined magnetic fields 
(dashed-black curves in Fig. [3E]). 

In addition, we have demonstrated (see Section 2|) 
that the <Syi? qu -criterion can recover, with much larger 
reliability than S/R quv , the original distributions of the 
three components of the magnetic field vector from the 
selected profiles. In particular, the latter criterion clearly 
overestimates the inclination (with respect to the observer's 
line-of-sight) of the magnetic field vector 7 (solid-red lines 
in Figs. HJi-HJl), and therefore also overestimates the com- 
ponent of the magnetic field that is perpendicular to the 
observer's line-of-sight, B± (solid-red lines in Figs. [TJj-Oj). 
To avoid these systematic errors, we propose employing 
instead the S'/i? qu -criterion, which allows the correct 
distribution of B± and 7 (solid-red lines in Figs.|3]3,d-|3jD,d) 
from the selected profiles to be retrieved. Unfortunately, 
as mentioned above, the S*/i? qu -criterion systematically 
selects the Stokes profiles where the magnetic field vector is 
already rather inclined, and thus ends up having the same 
sort of bias as the S/i? quv -criterion. It seems therefore 
that the only way around this problem is to decrease 
the photon noise to a sufficiently low level that the vast 
majority of profiles have a S/R ratio larger than 4.5 in 
the linear polarization profiles. In this way, the selected 
Stokes profiles would arise from magnetic fields that are 
representative of the real distribution. It is important to 
bear in mind that the overestimation of 7 and B± is a 
result that does not depend on the probability distribution 
function employed in our tests (Eq. [3]). 

We have also studied the ability of the inversion code 
to retrieve the angular distribution of the magnetic field 
vector (7 and c/>), even when the Stokes Q and U profiles 
are below the 4.5f7-level (Section [5]). In this case, we 
employed initial distributions of <f> and 7 featuring <5-Dirac 
functions. We found that the correlation introduced into 
the noise by these signals does not help us to retrieve 
the correct azimuthal angle <f>. Instead of the original 
(5-function, the inversion code yields a random distribution 
of azimuthal angles that are uniformly distributed (Fig. [5]). 
In the case of the inclination of the magnetic field vector 
7, the inversion code is able to partially retrieve the 
(5-function from the weak signals below the noise (Fig. [JJ). 
However, this is not the case when employing more general 
distributions, in which case the aforementioned systematic 
overestimation of 7, dominates the inferred distribution 
(cf. del Toro Iniesta et al. 2010). 

We caution that in our experiments we employed only 
one spectral line, Fe I (<? e ff = 2.5) 6302.5 A, whereas 
m many previous studies Fe I (g eS = 1.67) 6301.5 A 
had also been analyzed. Employing two spectral lines 
would certainly somewhat improve the results from these 
experiments as there would be more data points available 
to the inversion code. However, based on the results of 
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Orozco Suarez et al. (2010), these improvements will likely 
be of second order compared to the systematic errors 
presented in this paper. 

An additional point of consideration refers to the 
numerical code used to infer the magnetic field vector 
from the inversion of the radiative transfer equation for 
polarized light (VFISV). This code operates, with some 
modifications, under a Levenberg-Marquardt non-linear 
minimization algorithm (Borrero et al. 2010). We cannot 
exclude that other numerical schemes such as genetic 
(Lagg et al. 2004) or Bayesian algorithms (Asensio Ramos 
et al. 2008; Asensio Ramos 2009) are less affected by 
photon noise and therefore able to infer more accurately 
the distribution of the magnetic field vector. 
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Fig. 1. Results of the numerical experiment of the synthesis and inversion of Stokes profiles. The level noise added 
is a — 10~ 3 . Solid-black lines indicate the original distributions employed in the synthesis of 2 x 10 6 Stokes profiles 
(Equation [3]). Dashed lines show the distributions obtained employing only the Stokes profiles that are selected with the 
iS/ i? quv -criterion. The red lines display the distributions obtained from the inversion of the selected profiles: a) the absolute 
value of the component of the magnetic field vector that is parallel to the observer's line-of-sight |Bn| = |Bcos7|; b) the 
component of the magnetic field vector that is perpendicular to the observer's line-of-sight B± = B sin 7; c) the magnetic 
field strength or module of the magnetic field vector B; d) the inclination of the magnetic field vector with respect to 
the observer's line-of-sight 7; e) the azimuthal angle of the magnetic field vector on the plane that is perpendicular to 
the observer's line-of-sight cf>; f) the magnetic filling factor a. In all panels, the dashed-black and solid-red curves are 
normalized to the quotient of the number of profiles selected with each particular criteria to the total number of profiles 
synthesized (2 x 10 6 ). 
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Fig. 2. Same as Figure Q] but using a photon noise a = 3 x 10 4 and the 5 l /i? quv -criterion. 
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Fig. 3. Same as Figure Q] but using a photon noise a — 10 3 and the S , /i? qu -criterion 
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Fig. 4. Same as Figure Q] but using a photon noise a = 3 x 10 4 and the 5 l /i? qu -criterion. 
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Fig. 5. Mean values of the differences between the original physical parameters and the inferred values (AX^; circles) and 
standard deviations around the mean (dashed lines; cr x ): (panel-a) ABy and <TB n as a function of -Bp; (panel-b) AB± and 
ctb ± as a function of B±\ (panel-c) AB and <tb and a function of B; (panel-d) A7 and cr 7 as a function of 7; (panel-e) 
A(f> and as a function of 7; and finally, (panel-f) Aa and a a as a function of a. Blue shows the results obtained from 
a photon-noise level of cr = 10~ 3 and selection criterion S/R quv (from Fig. QJ; red represents a — 3 x 10~ 4 and S/R quv 
(from Fig. [2]); green is for a = 10 -3 and S/R qu (from Fig. [3]); and yellow is for cr = 3 x 10~ 4 and S/R qu (from Fig. [4}. 
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